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Abstract 

Using 13.53 ftr 1 of CLEO data, we have measured the ratios of the branching 
fractions fl+ = gfef^ , K = gg^ffi and the combined 
branching fraction ratio = B(D+^K-°Sn^) • We find i?+ = 0.74 ± 0.04 ± 
0.05, i?+ = 0.72 ± 0.10 ± 0.06 and Rf = 0.74 ± 0.04 ± 0.05, where the first 
and second errors are statistical and systematic, respectively. The known 
branching fraction B{D + -» R--k + t: + ) leads to: B{D + -» #*°e + z/ e ) = (6.7 ± 
0.4 ±0.5 ±0.4)%, £(L>+ -» i^*V + ^) = (6.5 ±0.9 ±0.5 ±0.4)% and 23(L>+ -> 
R*°l + vi) = (6.7 ± 0.4 ± 0.5 ± 0.4)%, where the third error is due to the 
uncertainty in B(D + — ► K~7r + 7r + ). 
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The transition amplitude of the semileptonic decay D + — > K*°l + vi is proportional to 
the product of the leptonic and hadronic currents. The hadronic current is represented 
by three analytic functions called form factors, Ax(q 2 ), A 2 (q 2 ), and V(q 2 ) JL|, where q 2 is 
the invariant mass squared of the virtual W boson. The form factors cannot be easily 
computed in quantum chromodynamics (QCD) since they are affected by significant non- 
perturbative contributions. Further, unlike the case of B — > D*lv, Heavy Quark Effective 
Theory (HQET) @-§J] cannot be applied in a straightforward way since the transition is to 
the light strange quark. Precise experimental measurements are needed to guide theoretical 
progress in this area. Measurements of these form factors are also valuable since, using 
HQET, they can be related to the form factors for the important decay b — > ulv, a heavy 
to light transition like the c — ► slv process we study here. This approach could reduce 
theoretical uncertainties on the value of \V u b\ measured in b — > ulv. 

In this work, the ratio of the branching fraction of D + — > K*°l + isi to that of D + — > 
K~7i + 7r + (Rf) was measured, where / stands for either an electron or a muon. In other 
experiments 0, the form factor ratios r 2 = ^ and ry = of D + — > K*°l + ui were obtained 
from angular correlations between the daughter particles. Combining these, we can calculate 
the form factors Ax, A 2 and V. 

The data used in this study was collected with the two configurations of the CLEO 
detector |J at the Cornell Electron Storage Ring (CESR). It consists of 9.13 fb _1 of integrated 
luminosity on the T(AS) resonance and 4.40 fb _1 below BB threshold. The investigated 
decay chain for this analysis was D* + — > D + n°, D + — > K*°l + i/, and K*° — > K~tt + . The 
use of D + produced in D* + — > D + ir° decays reduced the background. This analysis assumes 
that the thrust axis of the event approximates the direction of the D meson, as will 
be explained below. This approach does not work well for isotropic events. Thus, events 
whose ratio of the second and zeroth Fox- Wolfram moments || was less than or equal to 
0.2 were rejected. Specific ionization in the drift chamber and shower information in the 
electromagnetic calorimeter were used to select good electron candidates. Muon candidates 
were required to penetrate at least 5 interaction lengths in our muon detector. To select 
good 7r°'s, all 77 pairs for which |M 77 — M w o \ < 2.5a were accepted, where a is the standard 
deviation of the n° mass measurement, obtained as a function of 7i~° momentum using clean 
7T° data samples. Other kinematic criteria chosen to optimize our sensitivity are given in 
Table § 

To reconstruct the momentum of the v, two methods were used to obtain up to three 
values of p v . In the first method, one or two values for the v momentum were obtained 
assuming that the thrust direction of the event represents the D direction. Given pKni, 
the constraint that rriKniu = rriD+ provided an ellipsoid of allowed D momenta. We then 
took the two intersections between the ellipsoid and the direction of the D. When there 
was no intersection, the point on the ellipsoid that lay closest to the D direction was used. 
The values for the v momentum were then the difference between these D momenta and 
PkttI- The second method used the missing momentum of each event as an estimate of the v 
momentum. If the missing momentum gave a value of rriKniu much greater than the mass of 
the D + (we used mr)*+ for the limit for convenience), then the v momentum estimate from 
the second method was discarded. Among these three v momentum estimates, the one that 
gave the value of 5m = m Knlu7T o — rriKniu closest to the known value of m£>*+ — rn D + was 
chosen. According to our simulation, this selection method gave the v momentum closest to 
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TABLE I. Requirements on the kinematic variables. Here, 9 is the polar angle between the 
e + e~ axis and the momentum of the particle candidate. 



Variables 


Requirements 


\Ptt+\ 


> 0.5 GeV/c 


\pk-\ 


> 0.5 GeV/c 


\pA 


> 0.18 GeV/c 


\Pe\ 


> 0.7 GeV/c for |cos0| < 0.81 


Wn\ 


> 1.4 GeV/c for cos(9 < 0.61 




> 1.9 GeV/c for 0.61 < |cos0| < 0.81 


\PKitl\ 


> 2.0 GeV/c 




1.2-1.8 GeV/c 2 


\PKn\ 


> 0.7 GeV/c 



the generated v momentum 75% of the time. 

The previous CLEO analysis of this decay |L0| superseded by this study did not use 



v reconstruction. Instead it used the pseudo mass difference, 5m p = m K7T i n o — rriKni, and 
the resolution of 5m p was much inferior to 5m. The v reconstruction reduced the statistical 
uncertainty in the B(D + — > K*°l + v{) measurement by 20% based on studies using signal and 
generic continuum Monte Carlo (MC) events. However, since our method for choosing among 
the three solutions for \p v \ tends to bias the 5m distribution of the background toward the 
signal region, it is important to understand the contribution of the background that peaks 
in the signal region. The uncertainty coming from this effect was included in the systematic 
errors. 

Two quantities were fitted to obtain the number of signal events: m^n and 5m. First, the 
data were divided into 50 bins (25 bins for the smaller muon data) over the 5m range from 
0.135 GeV/c 2 to 0.235 GeV/c 2 . For each 5m bin, the Kn mass distribution was plotted and 
the number of K*'s was extracted by a fit. Second, these yields were plotted as a function 
of 5m, and the number of signal events was obtained by another fit. In the Kn mass fit, the 
signal shape was described by a p-wave K* Breit-Wigner distribution. The mass and width 
of K* obtained by the fit were consistent with the Particle Data Group (PDG) values ||. 
An analytic threshold function 

(m Kn - 0.64) a X e /3(^-0.64)+ 7 (m K . -0.64)' 

parameterized the shape of the background. The appropriateness of this function was tested 
using generic continuum MC events after discarding the signal among those events. The 
parameters a, j3 and 7 in this background function were unconstrained in the data fits. A 
typical Kti mass fit in the 5m signal region is shown in Fig. |I| 

The signal shapes for the 5m fits were obtained from signal MC samples, and the back- 
ground shape was described by a continuous function which was designed to accommodate 
the excess at m D *+ — m D + due to our neutrino momentum selection method. The shapes 
of these functions in the 5m fits for data were determined using a generic continuum MC 
sample. The 5m fit of the electron data is shown in Fig. Q. 
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FIG. 1. Kir mass distribution in the third 5m bin of the D + 
data. 
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6 



TABLE II. The electron, muon and electron-muon combined systematic errors in units of 10 



The combined systematic errors of the electron and muon modes are denoted as the I. 





Source 




e 


M 


I 


fits 
















5m sig. shape 


2.2 


3.3 


1.9 






6tti bgr. shape 




1 A 
1 .1 


i n 
1 . u 






data MC diff.(sig-) 


0.34 


3.8 


0.68 






data MC cliff. (bgr.) 


2.6 


2.5 


2.6 


lepton id 






1.1 


0.27 


0.89 


fake leptons 






0.064 


0.95 


0.17 


sig MC model dep. 






0.74 


0.72 


0.74 


efficiency by fragmnt. 






0.53 


0.54 


0.53 


feed down 






0.18 


0.17 


0.18 


D+ - K^^l+Vi 






2.6 


2.5 


2.6 


normalization mode 
















sig. shape 


1.1 


1.1 


1.1 






bgr. shape 


0.81 


0.79 


0.81 


TOTAL 






4.8 


6.0 


4.7 



The decay chain, D* + — > D + n°, D + — > K^n + n + , was used for the normalization mode. 
One advantage of this mode is that many systematic errors in its "measurement" are common 
with the D + — > K*°l + v mode and cancel in their ratio. Furthermore, the branching fraction 
of D + —>■ K~7t + tt + is well measured: (9.0 ± 0.6)% ||. Most of the kinematic requirements 
were the same as those in the semileptonic analysis, but the requirement on \pu* | was tuned 
to minimize our sensitivity to the uncertainty in the fragmentation of the charm quarks. 

For the normalization mode analysis, two successive fits of mr>+ and 5m were used. The 
data were divided into 50 5m bins from 0.135 GeV/c 2 to 0.185 GeV/c 2 . For each 5m bin, a 
D + mass fit was used to extract the number of D + 's. These yields of D +, s were plotted as 
a function of 5m. Then, this 5m plot was fitted to obtain the number of the normalization 
mode events. 

Table |TJ lists the important contributions to the systematic error. The dominant sources 
are the 5m signal shape estimation due to finite MC statistics, potential flaws in the estimate 
of the background shape obtained from the simulation of e + e~ — > qq and background arising 
from D+ -> i^VZ+z/j. 

The systematic error due to a potential difference in the background between the MC 
and data was estimated using a ir° mass sideband, 3.5a < \M~ n — M^o] < 15.0a. The ir° 
mass sideband data were analyzed as 7r° signal band data and the difference in yields of data 
and generic continuum MC sample was used in the estimation of the systematic error due 
to the deficiency of the background MC in reproducing the data. 

The upper limit of the branching fraction for D + — > K* q -k° ii + v^ with respect to D + — > 
K~ix + fJ^V/j, is 0.042 |ll| at the 90% confidence level. We used the product of this upper limit 
and the ratio of "efficiencies" that K*n Q lv and K*lv are detected as signal in our analysis 
as the systematic error due to this source. 

The systematic errors from each source were estimated individually for the electron and 
muon modes and then combined to obtain the overall lepton mode systematic errors. 

The ratio was calculated using 
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TABLE III. Results of the R£, Rt and R[ measurements and the branching fractions B e , B^ 
and Bi (see text). The first and second errors in the R and B measurements are due to statistical and 
systematic errors, respectively. The third errors in the B measurements are due to the uncertainty 
in B(D+ -> K-Tr + ir + ). 



R+ 


0.74 ±0.04 ±0.05 


K 


0.72 ±0.10 ±0.06 




0.74 ± 0.04 ± 0.05 


B e 


(6.7 ±0.4 ±0.5 ±0.4)% 


B fl 


(6.5 ±0.9 ±0.5 ±0.4)% 


B t 


(6.7 ±0.4 ±0.5 ±0.4)% 



TABLE IV. Comparison of the measured values of R£ and Ri 



Group 



R. 



Group 



CLEO(2£ 
PDG 

CLEO 

OMEG 
ARGTJSJp 

E691 O" 



0.74 ±0.04 ±0.05 
0.54 ± 0.05 
0.67 ±0.09 ±0.07 
0.62 ±0.15 ±0.09 
0.55 ±0.08 ±0.10 
0.49 ± 0.04 ± 0.05 



CLEO(2Q01) 
PDG M 
E687 2 
E653 [15 



0.72 ±0.10 ±0.06 
0.53 ± 0.06 
0.56 ±0.04 ±0.06 
0.46 ± 0.07 ± 0.08 



y-,-1- -^Vsi ^had /c-.\ 

where N 8 i and e s i are the observed number of events and the efficiency for D + — > K*°l + ui, 
respectively. Similarly, Nhad and etad are the corresponding quantities for D + — > K~tt + 7t + . 
The measured branching ratios Rt and Rf ', including their errors, are summarized in 



Table [II. The values of and R+ agree well. The PDG value of B(D + — > K 7r + 7r + ) = 
(9.0 ± 0.6)% H| was used to calculate the branching fractions of D + — > K*°e + v e (B e ), 
D + ^ R* ^ (Bp), and D+ - K*°l+v l {Bt). 

Table |TV] compares the results of this work with previous measurements [pUj-HTo]] and the 
PDG averages 0. Note that the E691 measurement of R£ and ours are the two most 
significant measurements, and they differ by about three standard deviations. We studied 
whether this difference (or part thereof) may arise from differences between the D + — > 
K*°e + v e model used here (ISGW2) @ and the one used in the E691 analysis (WSB) 0. 
Since the lepton momentum distribution expected from this decay has a significant effect 
on its detection efficiency in both experiments, we compared the electron energy and q 2 
distributions from the two models, but found no significant differences. We did not identify 
any other explanations for the discrepancy. 

Table |V| shows the decay rate for D + — > K* Q e + v e measured in this work and the the- 
oretical predictions p!6|JT8| -p2| . Most predictions are consistent with both our result and 
that of E691. Using our branching ratio measurements and the ratios between the form 



factors measured by E791 [0, we calculate Ai(0) = 0.69 ± 0.07, A 2 (0) = 0.48 ± 0.08, and 
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TABLE V. Comparison of the decay rate for D + — ► K*®e + v e measured in this work with those 
theoretically predicted and the present PDG value. All values are in units of 10 10 |I4 S | 2 sec -1 . 

Group Decay Rate (10 10 1 V ea | 2 s" 1 ) 

CLEO(2001) 6.67 ±0.75 

PDG § 4.79 ± 0.40 

UKQCD (2001) |l§ 5.8 ±0.5 

APE (l9) 6.9 ±1.8 

UKQCD (1994) 6.0+?;® 

ELC [EL 6.4 ± 2.8 

LMMS |gi 5.0 ± 0.9 

ISGW2 Ml 5.7 



V(0) = 1.25 ±0.15, where the errors are the quadratic sums of the statistical and systematic 
errors in this analysis and in the E791 measurement. 

The ratio of the vector and pseudo-scalar decay widths, r vp = rp^/re+iTf > tests quark 
models. This ratio was predicted to be in the range 0.9 to 1.2 by early quark models [17]|13 [27| 
and lattice gauge calculations f28| , |2"9|| . The E691 result implied that r vp = 0.48 ± 0.10. The 



newer ISGW2 model [IB| accommodated this measurement well, predicting r vp = 0.54. Our 
measurement of on the other hand, implies that r vp = 0.99 ± 0.06 ± 0.07 ± 0.06, where 
the first, second and third errors are due to our statistical and systematic errors and the 
uncertainty in B(D + — > K e + u e ), respectively, using the PDG value of B(D + — > K°e + i/ e ) = 
(6.8 ± 0.8)% 0. Our ratio is consistent with older models and disagrees with the ISGW2 
prediction by four standard deviations. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, the Natural Sciences and 
Engineering Research Council of Canada, and the Texas Advanced Research Program. 
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